Surface Criticality and Multifractality at Localization Transitions 



O 
O 

(N 

U 

Oh- 
< 

in 



A. R. Subramaniam 1 , I. A. Gruzberg 1 , A. W. W. Ludwig 2 , F. Evers 3 , A. Mildenberger 4 , and A. D. Mirlin 3 ' 5 '* 
James Franck Institute, University of Chicago, 5640 S. Ellis Avenue, Chicago, Illinois 60637, USA 
2 Physics Department, University of California, Santa Barbara, California 93106, USA 
3 Institut fur Nanotechnologie, Forschungszentrum Karlsruhe, 76021 Karlsruhe, Germany 
4 Fakultat fur Physik, Universitat Karlsruhe, 76128 Karlsruhe, Germany 
5 Institut fiir Theorie der kondensierten Materie, Universitat Karlsruhe, 76128 Karlsruhe, Germany 
(Dated: Received 2 December 2005; published 27 March 2006) 

We develop the concept of surface multifractality for localization-delocalization (LD) transitions 
in disordered electronic systems. We point out that the critical behavior of various observables 
related to wave functions near a boundary at a LD transition is different from that in the 
bulk. We illustrate this point with a calculation of boundary critical and multifractal behav- 
ior at the 2D spin quantum Hall transition and in a 2D metal at scales below the localization length. 
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Localization-delocalization (LD) transitions in disor- 
dered electronic systems represent a remarkable class of 
quantum phase transitions. One of their striking features 
is that moments of wave functions near a LD transition 
are characterized by a set of independent critical expo- 
nents 0] . Criticality of the wave functions manifests itself 
in their multifractality (MF) [2, 01 1 which has been stud- 
ied analytically and numerically for a variety of systems: 
Anderson transition ha d = 2 + e, 3, and 4 dimensions 
11 as well as weak MF random Dirac fermions 
(|, integer quantum Hall (IQH) 0, spin quantum Hall 
(SQH) 8] transitions in two dimensions, and power-law 
random banded matrices 

It is known that at conventional phase transitions, dif- 
ferent critical behavior occurs at a boundary, as com- 
pared to the bulk of the sample [lo|. This is especially 
well understood in two dimensions where methods of con- 
formal field theory (CFT) are available flll |. 

In this Letter we propose to extend the concept of 
surface criticality to LD transitions. Such an exten- 
sion has several potentially important applications. First 
of all, LD transitions are studied experimentally mainly 
through transport measurements, which are performed 
by attaching leads to the boundary of a finite sample. 
With such a setup it is feasible to measure surface criti- 
cal behavior directly. Secondly, we will show that even if 
the MF of wave functions in the whole sample is studied, 
the surface effects are fundamentally important. Further, 
LD transitions in two dimensions are expected to be de- 
scribed by certain CFTs. Such a description remains 
elusive for the primary example of the IQH transition, 
though some proposals have been put forward 0, 0] . 
Studying MF of wave functions near a surface is expected 
to help identify the CFT for the IQH transition and sim- 
ilar systems. 

Let us start with a brief review of bulk MF in the 
context of LD transitions. The MF of wave functions 
ip(r) at a LD transition is characterized by the scaling of 
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moments of ^(r)! 2 with system size L: 

L d mv)\ 2 «)~L~ T % r q = d{q- 1) + A q , (1) 

where (. . .) denotes the disorder average. Note than one 
often introduces fractal dimensions D q via,r q = D q (q—1). 
In a metal D q = d, while at a critical point D q is 
a non-trivial function of q, implying the MF of wave 
functions. Non-vanishing anomalous dimensions A q = 
(q — l)(D q — d) distinguish a critical point from a metal- 
lic phase and determine the scaling of wave function cor- 
relations. Among them, A2 < plays the most promi- 
nent role, governing the spatial correlations of the in- 
tensity |7/f: L 2d (|V 2 (r)V 2 (r')|) ~ (|r - r'\/L) A \ This 
equation, which in technical terms results from an op- 
erator product expansion of the field theory [J^, can 
be obtained from by using the fact that the wave 
function amplitudes become essentially uncorrelated at 
|r — r'| ~ L. Scaling behavior of higher order spatial cor- 
relations, (\tp 2qi (ri)i/) 292 (r 2 ) . . . ip 2qn (r n )\), can be found 
in a similar way. Above, the points were assumed to 
lie in the bulk of a critical system. In this case we denote 



the multifractal exponents by r q , A^, etc. 

In finite electronic systems wave functions vanish on 
the boundary. Therefore, to assign a meaning to the MF 
on the boundary, we have to interpret surface correlation 
functions in the sense of the theory of critical phenomena 
[ToL Hlj . Specifically, by "points on the boundary" we 
mean points that lie close to it, as compared to other 
distances, e.g. a few lattice spacings away. 

With this caveat in mind, we now generalize the notion 
of MF to systems with boundaries. First of all, we note 
that in general even the average value of IVK 1 ")! 2 f° r r 
at the boundary may be non-trivial, (^(r)! 2 ) oc L~ d ~ M . 
In particular, fi can be non-zero in systems with un- 
conventional (chiral or Bogoliubov-de Gennes) symmetry 
[lfij which are attracting particular interest in connection 
with physics of disordered superconductors. For the mo- 
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ments of the intensity, we get in analogy with Eq. , 

L*- 1 {\1>(t)\*>) ~ L-*l, (2) 
T s q =d(q-l) + qp+l + A s q , (3) 

with a new set of surface multifractal exponents t|, A^, 
which are in general independent on their bulk counter- 
parts 0| . The normalization factor L d ~ x is chosen such 
that Eq. (J2J) yields the contribution of the surface to the 
inverse participation ratios (P q ) — (J d d r\t/j(r)\ 2q ) con- 
ventionally studied in the framework of the MF analysis. 
The exponents A^ as defined in Eq. © vanish in a metal 
and govern statistical fluctuations of wave functions at 
the boundary, (|^(r)| 2 «)/{|^(r)| 2 ) 9 - L~ A «, as well as 
their spatial correlations, e.g. L 2 ( d+t ^{\i[) 2 (r)ip 2 (r')\) ~ 
(\r-r'\/L) A k 

In fact, these notions can be further generalized. First, 
for multi-point correlation functions, some points may be 
in the bulk while the rest on the surface. Then the scaling 
behavior will be described by mixed bulk-surface expo- 
nents. Second, surface and mixed exponents will sen- 
sitively depend on the global geometry of the boundary. 
For example, if the point r in Eq. J5J lies near the edge of 
a wedge with the opening angle 9, multifractal exponents 
will continuously depend on 9. We relegate the analysis 
of these generalizations to a future publication 0] and 
concentrate here on the fundamental surface exponents 
t s A s 

We now illustrate these points in the case of the 2D 
SQH plateau transition, which belongs to symmetry class 
C (in the classification of Ref. relevant for the de- 

scription of quasiparticle transport in singlet supercon- 
ductors with broken time-reversal symmetry A re- 
markable feature of the SQH effect is that a number of 
basic observables can be calculated exactly, as was dis- 
covered in Refs. [lflLlicl . via a mapping of the correspond- 
ing network model [21| (similar to the IQH-network of 
Ref. [2^]) to the problem of classical percolation. This 
mapping was extended in Ref. Q to extract analytical 
values of the bulk exponents A 2 = —1/4, A3 = —3/4. 

In what follows, we generalize this mapping to the SQH 
network with a boundary, and use it to extract new sur- 
face critical and multifractal exponents. The SQH net- 
work consists of directed links r that carry doublets of 
complex fluxes ip a (r) representing propagation of spin 
1/2 particles. Effects of disorder are introduced through 
random SU(2) scattering matrices on the links. At each 
node the scattering from two ingoing to two outgoing 
links is described by the (spin-independent) matrix 5, 
with S n = S22 = (1- i 2 ) 1/2 , £12 = -Sai = t. The value 
t = l/\/2 corresponds to the critical point of the SQH 
transition. 

In order to study the effect of boundaries, we impose 
reflecting boundary conditions along one direction. This 
does not affect the symmetry of the system, which en- 
ables us to retain the mapping to percolation. In tech- 



nical terms, within the approach of Ref. |19j , supersym- 
metry (SUSY) is preserved at the boundary nodes of the 
network model. We have also checked that the SUSY 
method can be extended to the two- and three-point func- 
tions needed to extract r q with q = 2,3, and yields re- 
sults identical with the path integral approach developed 
in Refs. HII|. 

We calculate first the average local density of states 
(LDOS) at a point ri on the boundary that will allow us 
to find the average of the intensity |?/> a (ri)| 2 . The LDOS 
can be expressed in terms of one-point Green's functions 
and becomes, when mapped to percolation, 

(p(r 1 ,e)) = (l/27r)[l-5^P(ri;JV)co8 2JVe], (4) 

JV 

where e is the energy and P(ri; N) is the probability of a 
iV-link hull passing through ri , in analogy with [2(J ■ The 
corresponding surface critical exponent is xf = 1/3 [Ijjj], 
which should be contrasted to its bulk value x\ = 1/4. 
The latter value implies that the percolation hull has 
fractal dimension 2 - x\ = 7/4, so that P(r, N) ~ A" 8 / 7 
for r in the bulk. This yields, according to Eq. (J2J, the 
DOS scaling p(e) cx e^/^i) = e 1 / 7 UlEl. Note that 
2 — x\ = 7/4 is the dynamic exponent governing the scal- 
ing of energy with the system size L at SQH criticality, 
so that the level spacing at e = (and thus the charac- 
teristic energy of critical states) is S ~ L~ 7 / 4 . 

In our case, when the point ri is located at the surface, 
we find P(t u N) ~ N~ 1 ~ x '^ 2 ~ a ^ = A~ 25 / 21 and the 
LDOS scaling p(ri,e) cx e x i/( 2 ~ x i*> = e 4 / 21 . For the wave 
function at the boundary, this implies 

L 2 (|Vari)| 2 )~£^=£- 1/12 . (5) 

Therefore, the average intensity of a critical wave func- 
tion is suppressed at the boundary with the exponent 
/i = 1/12. A similar calculation for the conductance be- 
tween two point contacts ri and r 2 located at the bound- 
ary [24| yields scaling with exponent 2x\ = 2/3, 

^ t (r 1 ,r 2 )oc|r 1 -r 2 |- 2 / 3 . (6) 

We turn now to the multifractal exponents. To calcu- 
late A|, we consider the correlation function ||: 

(2^)- 2 P(r 1 ,r 2 , ei , £2 ) = (]T llMrOfhMra)! 2 

x 6(e! - 6i)S(e 2 - ej)}. (7) 

To study critical states, we take ei j2 = and broaden 
delta-functions by 6. Using the mapping to percolation, 
we relate (|7J| to percolation probabilities, 

P(n, r 2 , z) = (1 - z 2Ar )P(i"i, r 2 ; N) 

N 

+ 4^ (l-z™)(l-z™')P-( ri ,T 2 ;N,N'), (8) 

N,N> 
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in analogy with Ref. j8|. Here z 



P(n,r 2 ;JV) is 



the probability of an iV-link hull passing through ri,r2, 
and P_(i"i, r 2 ; N, N') is the probability of an TV-link hull 
passing through ri and a different iV'-link hull passing 
through r 2 . In view of the cancellation of leading order 
terms at z — ► 1, we need to consider the next, sub-leading 
term. Using again the surface critical exponent x\ = 1/3, 
we find, when both ri and r 2 lie at the boundary 



-25/21-1/3 



ri - r 2 



Substituted in Eq. (JSJ, this gives the value 
A* = -x\ = -1/3. 



I < N^ 7 . 



(9) 



A similar calculation can be done for the three-point 
function, with the result 



A s — 

^3 — 



-3x1 



-1. 



(10) 



In this case there is no cancellation in the limit z — > 1, 
and the result conforms with the usual scaling of a three- 
point function at the boundary. 

We have verified these predictions by direct simulation 
of the SQH network. To this end we have numerically 
diagonalized [25| the AL 2 x AL 2 discrete time-evolution 
operator of the network model and selected, for each re- 
alization of disorder, two eigenstates ipi a (r) with low- 
est quasi- energies e*. We found from the scaling of 
the moments (IV^q ( r )\) averaged over the ensemble of 
10 4 — 10 5 systems. As shown in Fig. ^ the numerical 
data fully confirm the values of the exponents fi, A|, 
and A| obtained above, Eqs. ©, ©, and JUJ- Clearly, 
numerical simulations allow us to study also the scaling 
away from the analytically accessible values q = 2, 3. In 
the lower panel, we show the results for the exponents 
A^ [divided by q(l — q)] in the range of q between and 
3.5 and demonstrate their independence from the bulk 
exponents A£. 

We also show in Fig. [3 the singularity spectrum /(a), 
obtained from r(q) by Legendre transformation. Its 
meaning is as follows: the number of points r in the sam- 
ple where the wave function intensity is |^ 2 ( r )l ~ L~ a 
scales as L^ a \ The difference between the maximal val- 
ues (2 vs. 1) of / b (a) and f s (a) simply reflects the dif- 
ferent dimensionalities of the bulk and the surface. On 
the other hand, the difference in the position of the max- 
imum («q ~ 2.137 vs. «g ~ 2.326) and in the width of 
the curve demonstrates that at the boundary the typical 
value of the intensity is suppressed, while fluctuations are 
stronger than in the bulk. 

We have thus shown that surface MF differs signifi- 
cantly from that in the bulk. One can now ask the follow- 
ing question. Imagine that one performs a multifractal 
analysis for the whole sample, without separating it into 
"bulk" and "surface" . Would then the surface exponents 
play any role? A naive answer is no: since the weight of 
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FIG. 1: Surface MF at the SQH transition. The numeri- 
cal data confirm analytically found exponents (dashed lines), 
with deviations as low as 4% (/i), 0.5% (A|) and 6% (A|). 
Lower panel: numerical results for the surface MF spectrum 
Aq for < q < 3.5. For comparison, the bulk spectrum A, is 
also shown, as well as the analytical results for q — 2, 3. 

surface points is down by a factor 1/L, one could expect 
that only the bulk exponents would matter. This is not 
true, however. To illustrate this point, we turn to an 
example, where the whole multifractal spectrum can be 
studied analytically. Specifically, we will consider a 2D 
weakly localized metallic system (dimensionless conduc- 
tance g^l), which shows weak MF [j| on length scales 
below the localization length £ ~ e (™9) j where (3 = 1 (2) 
for systems with preserved (broken) time-reversal sym- 
metry. With minor modifications, the formulas below 
describe also the Anderson transition in 2 + e dimensions. 

The bulk multifractal spectrum of this system was ob- 
tained via the perturbative renormalization group treat- 
ment of the underlying field theory (er-model) and also 
within the instanton approach [5j. The result reads 

T h q =2{q-l) + 1 q{l-q)- J = (fag)' 1 « 1. (11) 

Generalizing this calculation to the surface case, we find 



T s q =2(q-l) + l + 2 lq (l-q). 



(12) 



The factor of 2 in front of the last (anomalous) term 
can be traced back to the corresponding enhancement of 
the return probability near the surface. Performing the 




FIG. 2: Surface and bulk f(a) spectra at the SQH transition 
obtained from the data of Fig. Q 
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FIG. 3: Surface and bulk multifractal spectra r q and f(a) for 
a 2D metal with 7 = 0.01. For details see text. 

Legendre transformation, we find the /(a)-spectra, 

f b (a) = 2-(a-2- 7 ) 2 /4 7 , (13) 
/» = l-(a- 2 - 2 7 ) 2 /8 7 . (14) 

These results are illustrated in Fig. [21 When the MF 
in the whole sample is analyzed, the lowest of the T q 
exponents "wins" . It is easy to see that the surface effects 
become dominant outside the range q_ < q < q + , where 
q± ~ ± 7 -1 / 2 are the roots of the equation = t|. The 
lower panel of Fig.[3]shows how this is translated into the 
f(a) representation. The total singularity spectrum is 
given by the bulk function / b (a) only for a b < a < a_, 
where a± — 2 ~ =f2 7 1 / 2 . Outside this range the surface 
effects are important. Specifically, f(a) is equal to the 
surface spectrum f s (a) for a < a s + and a > evl, where 
a± — 2 ~ =f4 7 1 / 2 , while in the intermediate intervals 
a s + < a < a\ and a^L < a < a s _ its dependence on 
a becomes linear (shown by dashed lines). The latter 
behavior is governed by intermediate (between "bulk" 
and "surface") points with a distance from the surface 
r r~i L 13 , < (3 < 1; their f(a) spectrum is easily found 
to be f p (a) = fif h (a) + (1 - /3)/ s (a). Note that in this 
case the surface effects modify f(a) in the whole range 
below /(a) ~ 1. Therefore, the surface exponents affect 
the multifractal spectrum of the sample not only for rare 
realizations of disorder (governing the negative part of 
f{a)) but also in a typical sample. 

To summarize, we have developed the concept of sur- 
face MF for localization transitions in disordered elec- 
tronic systems, and have extended the notion of surface 
criticality to those transitions. We have calculated sur- 
face critical and multifractal exponents for the SQH tran- 
sition. Considering the example of a 2D weakly local- 
ized system, we have shown that the surface criticality 
may crucially affect multifractal spectra in systems with 
boundaries. Our work opens a new direction of research 
in the field of Anderson and quantum Hall transitions. 
More generally, it is interesting to study boundary ef- 
fects for MF in other stochastic systems. 
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